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(I)  INTRODUCTION 


Transonic  channel  flows  have  applications  in  most  modern  high 
performance  jet  engines.  Because  the  blades  in  a rotor  rotate  past 
fixed  vanes,  and  because  of  various  disturbances  caused  by  combustion 
instabilities,  flow  asymmetries,  etc.,  these  flows  are  inherently  un- 
steady. Since  the  flow  is  transonic,  small  changes  in  pressure  down- 
stream of  a shock  wave,  can  cause  large  variations  in  the  position  of 
the  shock  wave  and  thus  large  variations  in  forces  on  the  blades. 

Previous  work  on  unsteady  two  dimensional  channel  flows  with 
shock  waves  indicates  that  a great  many  physical  problems  of 

interest  fall  in  the  so  called  slowly  varying  time  regime.  If  is 

the  characteristic  time  associated  with  an  impressed  pressure  oscil- 
lation,  say,  (overbars  denote  dimensional  quantities)  and  L and  a are 
the  channel  throat  half  width  and  critical  sound  speed  respectively, 
then  this  region  is  characterized  by  the  relation 


a 

That  is,  the  period  associated  with  the  impressed  oscillation  is  large 
compared  to  the  time  taken  by  an  acoustic  wave  to  cross  the  channel. 
As  pointed  out  in  reference  I,  for  typical  conditions  in  a jet  engine. 


'^^Richey,  G.  K.  and  Adamson  Jr.  , T.  C.  . "Analysis  of  Unsteady 
Transonic  Chaiinel  Flow  with  Shock  Waves,  " AIAA  .Tournal,  14,  197(>, 
pp.  1054-1061. 

^^^Messiter,  A.  F.  and  Adamson  Jr.  , T.  C.  , "Asymptotic  Solutions 
for  Nonsteady  Transonic  Channel  Flows,  "Symposium  Transsonicvim  II, 
Eds.,  Oswatitsch,  K.  and  Rues,  D.  , Springe  r Verlag,  1976,  pp.  41-48. 
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this  "slowly  varying"  time  regime  encompasses  a range  of  frequencies 

(f  ■ ~ cycles/sec)  of  100  < f < 1000.  This  frequency  range 

is  evidently  of  considerable  interest. 

Within  the  general  condition  T . /(L/a  ) » 1 there  are  two 

ch 

specific  i;ases  of  physical  interest  in  flows  with  shock  waves.  They 

have  to  do  with  the  comparison  of  the  characteristic  time  associated 

with  the  impressed  pressure  oscillations,  T , , with  the  characteristic 

ch 

time  for  a signal  to  travel  from  the  point  where  the  oscillation  is  intro- 
duced downstream  of  the  shock  wave  to  the  shock  wave  itself,  say 

T , Thus,  if 
sh 

^ch 

- 0(1)  pressure  signals  from  downstream 

sh 

reach  the  shock  wave  in  a time  com- 
parable with  the  time  characterizing 
the  pressure  variation  e.g.,  the 
period  of  the  pressure  oscillations  ). 
Oscillations  in  velocity,  pressure, 
etc.  , in  the  channel  flow  downstream 
of  a shock  wave  lag  behind  the  im- 
pressed pressure  oscillations. 


(2) 


» 1 


pressure  signals  from  downstream 
reach  the  shock  wave  "inst antaneous - 
ly"  compared  to  the  time  character- 
izing the  pressure  variation.  Oscil- 
lation.s  in  velocity,  pressure,  etc., 


L 


in  the  channel  flow  downstream  of  a 


shock  wave  are  in  phase  with  the 
pressure  oscillations. 

The  remaining  order,  « 1.  appears  to  be  a more  complica- 

ted and  probably  less  important  version  of  case  (1),  and  is  not  con- 
sidered here.  Case  (1)  is  considered  in  reference  (1)  for  a symme- 
tric channel.  In  reference  (2),  solutions  are  presented  for  case  (1) 
written  in  terms  of  unsteady  boundary  conditions  rather  than  initial 
conditions;  in  addition,  the  concepts  inherent  in  case  (2)  are  dis- 
cussed briefly  and  the  differential  equation  describing  the  instantane- 
ous shock  wave  position  in  a channel  with  parabolic  walls  is  derived 
for  a particular  case  in  which  the  impressed  disturbances  have  a 
sinusoidal  time  dependence.  In  both  references,  solutions  are  given 
in  term.s  of  asymptotic  expansions;  thus,  the  type  of  solution  first 
postulated  by  Szaniawski  for  steady  transonic  flows  are  derived 
in  a systematic  manner  and  extended  to  cover  unsteady  flows. 

In  this  report,  solutions  for  case  (1)  are  reviewed  and  numeri- 
cal solutions  for  a typical  case  are  given.  Next,  detailed  solutions 
for  case  (2)  are  presented,  again  with  numerical  results  for  a typical 
case.  Finally,  a general  solution,  (case  (3)),  which  holds  for  either 
case,  is  presented,  with  an  error  estimate  which  allows  one  to  make 

^Szaniawski,  A.,  'Transonic  Approximations  to  the  Flow  Through 
a Nozzle,  " Archiwum  Mechaniki  Stosowanej,  17,  1965,  pp.  79-85. 
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simplifications  in  the  numerical  calculations  on  a systematic  basis, 
i.e.  , on  the  basis  of  a desired  accuracy.  A detailed  description  of 
the  various  possible  shock  wave  motions  for  periodic  impressed  pres- 
sure oscillations  is  also  presented.  All  results  are  given  for  symme- 
tric channels. 

Finally,  a discussion  of  the  use  of  steady  flow  predictions  of 
shock  induced  separation  to  predict  separation  in  unsteady  flows  is 
given. 


(II)  UNSTEADY  TRANSONIC  CHANNEL  FLOWS 
WITH  SHOCK  WAVES 

The  problem  considered  is  that  of  a two  dimensional  channel 
with  an  arbitrary  wall  shape,  in  which  there  is  a transonic  flow  with 
a shock  wave.  A sketch  of  the  channel  flow  illustrating  the  notation 
used,  is  shown  in  figure  1.  The  flow  upstream  of  the  wave  is  steadv^; 
pressure  oscillations  are  impressed  upon  the  flow  downstream  of  the 
shock,  at  X = X,  say,  causing  the  shock  wave  position  to  oscillate. 

The  gas  is  assumed  to  follow  the  perfect  gas  law  and  to  have  constant 
specific  heats.  The  flow  upstream  of  the  shock  wave  is  isentropic, 
and  because  the  flow  is  transonic,  the  shock  is  weak  enough  that  a 
velocity  potential  may  be  used  to  the  order  desired.  The  coordinates, 
X and  y,  are  made  dimensionless  with  respect  to  L,  the  time,  T,  with 
respect  to  L/a  , and  velocities  with  respect  to  a ; hence,  the  di- 
mensionless velocity  potential  is  referred  to  the  product  La  . The 
pressure,  P,  density,  p,  and  temperature,  T,  are  made  dimension- 
less with  respect  to  their  critical  values,  and  the  enthalpy  is  referr(;d 


4 


p 


-*  2 

to  a 

The  wall  shapes  considered  are  written  as  follows,  for  symme- 
tric channels : 

= + (1  + f(x))  (1) 

where  f(x)  is  the  arbitrary  wall  shape  function,  such  that  f(0)  = f'(0) 

= 0.  Thus,  X is  measured  from  the  channel  throat.  The  radius  of 

_2 

curvature  of  the  channel,  at  the  throat,  is  0(r  ) from  eqn.  (1),  and 

as  will  be  seen  later,  u-i  = 0(e).  Hence,  for  transonic  flow,  e « 1. 

Unsteady  flows  may  be  characterized  by  prescribing  the  rela- 
tive order  of  the  characteristic  time  associated  with  the  impressed 

disturbances,  T , and  the  characteristic  time  associated  with  the 
c n 

_ 5j: 

acoustic  waves  traveling  through  the  channel,  L/a  . As  mentioned 
previously,  the  so  called  slowly  varying  time  regime  is  considered 
here,  where  T^  » L/ a . Therefore,  a parameter,  t , is  intro- 
duced and  the  time  is  stretched,  as  follows: 

T 5 t^j^/(L/a')  T ^ T t (2a,  b) 

so  that  T » 1 and  t = 0(1).  The  relationship  between  t and  c depends 
on  which  of  the  cases  mentioned  previously  is  considered.  Again, 

T^h  order  of  the  time  it  takes  a signal  to  travel  upstream  from 

the  origin  of  the  flow  disturbance  to  the  shock  wave,  a distance  of 
order  L,  say.  A disturbance  pulse  travels  upstream  at  sonic  veloc- 
ity relative  to  the  flow  and  thus  at  an  absolute  dimensionless  velocity, 

u„  “ 1 -u  = Ole).  Hence,  T , - 0(L/a  u_,)  = 0(L/a  t),  and 

sh  P 

= 0(t/c  ).  Therefore,  it  is  seen  that  case  (1)  mentioned 
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previously  and  considered  in  reference  (1)  is  that  for  which  t-  0(c  ), 

- 1 -2. 

and  case  (i)  is  that  for  which  t » c , and  here  we  choose  t - 0(r  ). 

In  summary,  then,  we  consider  the  following  two  cases. 


Case 

(1) 

t3al 

Case 

(2) 

(3b) 

where  k,  and  k;>  are  arbitrary  constants  of  order  unity. 

1 ^ 

In  the  following,  solutions  are  presented  in  the  form  of  asymptotic 

2 

expansions,  uniformly  valid  to  order  c . Although  the  method  is  valid 
for  more  general  impressed  disturbances,  attention  here  is  focused  on 
oscillatory  second  order  pressure  fluctuations  impressed  on  the  flow 
downstream  of  the  shock  wave.  Upstream  of  the  shock  wave,  the  flow 
is  steady. 


Solutions  for  Case  (1),  x ^ = (k^  e) 


-1 


The  solutions  for  Case  (1),  derived  in  reference  (1),  may  be  writ- 
ten as  follows  : 


1+cUj+e  u^+... 

(4a) 

2 

c v^  . . . . 

(4b) 

2 

1 - eyuj  - r • 

(4c ) 

1 - CUj  - ^1' 

) + . . . 

(4d) 

2 

1 - c(y-  1 )Uj  - 1 Ku^  + 

(4e ) 

where  v = C,.,/C  is  the  ratio  of  specifii'  heats,  and  where 
• I P V 


- f(x)  t C 
1 ) w 


( *ia  I 


b 


(5b) 


X 


V = f'y  + C 


(5c) 


4f';l(Y+l)C  ™ n ..  , 

j — Z/  — ^ expf-nTTx  /1(y+1)C^]  } 

n n = l n 


cos  (n  7T  y ) 


(x-x  )e 
o 


-1/2 


(5d) 
I F>e ) 


i If'  (-  (2Y-3)uf]  + -^  + -Z  G(t-t  ) 
6 I u j f 


(5f) 


~ * J (2kj /(y^I  )Uj  ))d| 

X 


(5gl 


Finally,  the  shock  wave  location  is 

x=x„+ex,(t)-t...  (b) 

s 0 si 

where  x_  is  the  steady  state  location,  and  x is  found  from  the  equa- 
0 s 1 

tion. 


4kj  dxs j 

f " 

_ h fh 

(7a) 

(Y+1)  dt 

X , X 3 u 

d u 

t - T / T j 

(7b) 

The  y dependence  of  x^  occurs  first  in  a term  of  order  e In 

eqn.  (5a),  is  an  arbitrary  constant  set  by  the  value  of  the  velocity 

at  the  throat;  i.e,,  if  the  flow  is  supersonic  or  subsonic  Ihei  e, 

C >0,  while  if  the  flow  is  one  which  accelerates  from  subsonic  to 
w 

supersonic  flow,  C =0.  Also,  in  eqns.  (5),  f df/dx,  etc.  , and  f. 
^ w ’ 

is  the  derivative  of  f evaluated  at  Xq,  is  the  value  ol  Uj  at  x^ 


t'valuatfd  upstream  of  tlie  shock  (upper  sign  in  eqn.  (5a)),  is  an 
arbitrary  c onstant  of  inU'gration  set  by  Ijoundary  conditions,  and  G(t) 
is  the  oscillation  in  pressure  impressed  downstream  of  the  shock  wave 
at  X - X,  say.  Hence,  upstream  of  the  shock  wave,  the  equation  for  h^ 
w(3uld  be  found  from  eqn.  (5f)  l)y  setting  G = 0.  In  cqns.  (5f)  and  (5g), 
t^  - ' |,(x)  is  the  lag  time,  representing  the  time  taken  by  a signal  to 
travel  upstream  from  x X.  to  any  position  x.  The  sul'scripts  u and  d 
in  eqn.  |7a)  indicate  functions  evaluated  immediately  upstream  and  down- 
stream of  the  shock  wave,  respectively;  thus,  for  example,  = u^^ 


and  from  the  shock  jump  conditions 


"id  = *^u 


(8) 


The  function  <,  '(x',y)  ineqns.(5)  is  the  additional  velocity  poten- 

I 

tial  needed  to  make  the  solutions  uniformly  valid  to  order  e . That  is, 
the  solutions  with  ^ - 0 are  the  outer  solutions,  which  do  not  satisfy 

the  shock  wave  jump  conditions  to  second  order.  Hence,  an  inner  re- 
gion enclosing  the  shock  wave  must  be  i onsidered.  In  this  inner  region, 

1/2 

which  has  an  extent  of  order  c so  that  the  inner  x-direction  coordi- 

nate, x^,  is  written  as  in  eqn,  (5e),  it  is  found  that  the  outer  solutions 
do  hold  upstream  of  the  shock  wave,  but  that  a new  solution  is  neces- 
sary downstream  of  the  shock  wave.  This  new  inner  solution  satisfies 
the  shock  junip  conditions  at  x - 0 and  matches  with  the  downstream 

outer  solution  as  x — m.  The  inner  and  outer  solutions  are  then  joined 

2 

to  form  a composite  solution  uniformly  valid  to  order  c . 1 hus  , up- 

stream  of  the  shock  wave  (,  ~ K - 0 in  eqns.  (5b)  and  (5c),  while  t 

and  C sire  found  from  eqn.  t5d)  downstream  of  the  sliock. 

y 


H 


Numerical  solutions  are  carried  out  for  (jivt'n  functional  forms 


for  the  wall  shape,  f(xl,  and  impressed  pressure  distribution,  G(tl,  and 

for  given  values  of  y,  c,  t,  and  the  steady  state  pressiu'e  at  x X. 

Again,  x = X is  the  point  at  which  a given  pressure  is  impressed  upon 

the  flow;  setting  a given  steady  state  pressure  there  sets  x^,  the  steady 

state  shock  le.c  ation,  and  prescribing  G't)  there  allows  calculation  of  the 

corresponding  nonstationary  shock  wave  location.  It  should  be  noted 

2 

that  variations  of  pressure  of  order  r at  x X,  are  sufficient  to  cause 
variations  in  shock  location  and  thus  pressure  variations  of  order  < in 
the  region  in  which  the  shock  wave  travels.  This  point  is  illustrated  in 
figure  2,  where  a typical  pressure  distribution  for  an  accelerating  flow 
has  been  sketched. 

Example  calculations  are  shown  in  figures  3 and  4 for  the  follow- 
ing conditions : 


> 


f(x) 

= 18x"/13 

X 

< 1 

T 27tx-2)‘^/)3 

» 48(x-2)'Vl3  + 3 

1 < X 

< 2 

= 3 

X 

> 2 

(9) 

G(t) 

- 3 s in  ( 2 f ) 

t 

> 0 

- 0 

t 

< 0 

r = 

0.  I 

Y ^1.4 

^2d 

- - 

2x  C V3 
' u 

T 

20  C = 

w 

0 

X = 3 

^2u 

0 

The  flow  is  taken  to  be  an  accelerating  flow  which  goes  through  sonic 

velocity  at  the  throat.  Hence  C =0  and  C,  ~ 0.  For  C 0,  satis - 
' w 2u  2u 

faction  of  the  shock  jump  conditions  demands  that  - ZyC^^/i; 


9 


I 


giving  a numerical  value  for  then,  is  equivalent  to  setting  the 

pressure  at  x = X and  sets  x^,  the  shock  location.  That  is,  if 


ix 


= - 4.  81^3 


tlO) 


then,  using  the  equations  given  for  f(x)  in  equations  (9),  one  can  show 
that  Xq  = 1.5.  Finally,  for  the  given  values  for  t and  c,  it  is  seen 
from  eqn.  (3a)  that  = 2.  With  these  parameters  known,  eqns.  14)  - 
(7)  may  be  used  to  calculate  velocity,  pressure,  temperature,  and 
density  distributions,  and  shock  wave  location  and  velocity,  where 


s ^ SI 

= k . f -T t 


(11) 


s dT  ' r dt 
It  should  be  noted  that  the  wall  shape  shown  in  eqn.  (9).  was  i hosen  so 
that  ftx),  f'(x),  and  f”(x)  are  continuous  throughout  the  channel. 

Isotachs  and  pressure  distributions  at  the  centerline  and  at  the 
wall,  at  given  times,  are  shown  in  figures  3a  - 3f  for  the  case  where 
the  flow  is  steady  until,  at  time  t = 0,  the  pressure  oscillation  indicated 
in  eqn,  (9)  is  started  at  x = X.  It  is  seen  that  as  the  disturbance  propa- 
gates upstream,  the  velocity  profiles  downstream  of  the  shock  change, 
and  at  time  t = 0.39456  the  disturbance  reaches  the  shock  and  the  shock 
position  begins  to  change.  The  farthest  forward  position  of  the  shock 
wave  occurs  at  t = 1.96536,  figure  3f. 

Centerline  pressiiiu;  distributions  corresponding  to  three  dif- 
ferent shock  locations  and  thus  to  several  different  times  are  shown 
in  figure  4.  The  shock  is  located  at  its  farthest  downstream  location, 

^ - 1.5,  until  t - 0.  39456;  thus  only  the  pressure  dislriliution 


1 0 


X = X 


downstream  of  x - 1.5  varies  with  time  as  illustrated  by  the  two  dis- 
tributions at  t 0 and  t - 0.  39456.  As  the  shock  bcj^ins  to  move  up- 
stream. the  fluid  velocity  relative  to  the  wave,  u - u , increases 
’ ' us 

because  the  increase  in  -u  , the  shock  wave  velocity  in  the  direction 

s 

of  the  throat,  is  greater  than  the  decrease  in  u^,  the  steady  state 

velocity  upstream  of  the  shock  wave,  caused  by  the  decrease  in  the 

cross-sectional  area  of  the  channel.  Hence  the  pressure  jump  across 

the  shock  increases,  as  indicated  by  the  distribution  marked  t - 1.  17996, 

where  x - 1.  18.  At  the  farthest  upstream  shock  location,  indicated 
s 

by  the  distribution  market  t = 1.  96536,  the  shock  velocity  is  zero. 

Finally,  as  the  shock  moves  downstream  and  again  reaches 

X = 1.38,  at  t = Z.  75075,  the  fluid  velocity  relative  to  the  shock  wave 

s 

is  less  than  the  steady  state  velocity  at  that  point.  Hence,  the  pres- 
sure jump  across  the  shock  and  indeed  the  pressure  at  any  point  dcwn- 
stream  of  the  wave,  is  less  than  the  pressure  at  the  corresponding 
point  for  the  I - 1.  17997  distribution  when  the  wave  was  advancing. 

In  figure  4,  it  should  be  noted  that  the  fact  that  several  of  the  pressure 
distributions  have  common  points  at  x = 3,  is  fortuitous. 

In  figure  5,  the  shock  wave  position  and  velocity  are  shown  as 
functions  of  time,  for  this  sample  case. 

Solutions  for  Case  (Z),  ^ ' 

The  fundamental  ideas  underlying  the  method  of  solution  for  this 
case  were  discussed  briefly  in  reference  (Z).  Ih're,  detailed  solutions 
are  given  for  the  same  general  problem  considered  foi  the  case  (1  I 

1 1 


ah'ulations . That  is,  unsteady  flow  with  shock  waves,  in  a two  dimen- 


sional channel,  is  considered.  The  unsteadiness  arises  as  the  result 
of  oscillations  impressed  on  the  flow  at  some  point  downstream  of  the 
shoi'k;  upstream  of  the  shock,  the  flow  is  steady.  '1  he  impressed  os- 
1,' illations  again  have  an  amplitude  of  order  < , i.e.  they  occur  in  the 


second  order  pressure.  Now,  however,  the  period  of  the  oscillation 
is  large  c ompared  to  the  time  taken  for  a .signal  to  travel  upstream 
'1  h . k wave . 

Bei  ause  ^ t ’ partial  derivative  with  respect  to 

time  is  ordered  by  t * fi.e.,  9 /9  T = ^ it  is  clear  that  insofar 

as  the  general  form  of  the  outer  solutions  for  u,  v,  P,  etc.  , are  con- 
cerned, the  solutions  for  case  (2)  should  be  derivable  from  those  for 


case  1),  i)y  simplification.  Thus,  eqns.  M ),  (5a),  and  eqns.  (5b)  and 
(5c)  w'ith  c,  0,  hold  for  tliis  case  also.  7 )ie  questions  which  remain, 
tlien,  are  tliosi-  concerning  the  calculations  in  the  inner  region,  for  < , 
tile  calculations  for  h , and  the  calculations  for  the  shot,  k wave  position. 

In  case  (1),  because  unsteady  oscillations  in  pressure  (and  thus  u) 
are  impresstul  (jnly  in  second  order,  the  first  order  solution  downstream 
of  l)ie  shock  has  to  be  (tie  steady  sulisonic  solution  asso< dated  with  sonic- 


flow  at  tlu-  throat,  i.e.,  u,  ~ - 


f'xl  (C  0),  since'  it  would 

( yt  1 ) w 


.1  liidd  for  the  case  wtiere  no  shock  e\ists,  as  indicated  in  figure 
;.  Itetwc'cn  the  throat  and  the'  slirick,  on  the  other  liand,  u^  '1  ^ 


I lence  u , 


u . As  a result  of  these  considerations,  the  shock 
I u 


'locity  must  be  of  order  c I hat  is,  the  sliock  wave'  jetnp  condition 


'lativc  to  file  wave  is,  to  first  order. 


12 


Id  Is 


(u  - u,  ) 

1 u Is 


where'll  fs  the  iirst  orde r-s+it)ck  velocity.  Since  u,  ,=  -u,  , it  is 

Is  ^ Id  lo 

seen  that  u =0,  i.  e.  the  shock  velocity  is  of  second  order.  Since 
I s 

Tj  0(c  S.  th(>n,  thi-  shock  velocity,  u^,  is 
dX  dX 

",  ' •'I'dT  = S " 5T  '■‘o  """i 

so  that  X must  be  constant,  and  x , = x ,(t).  In  the  present  case, 
o s 1 s 1 

for  the  reasons  mentioned  above,  the  shock  wave  velocity  is,  again, 

I ^ -2 

of  order  r . However,  = 0(c  ),  so 

dX  , dX 

s , ^ s , , , , 

"s  = ^ 

and  we  see  that  now  x varies  with  t in  zero  ^ order.  That  is,  the 

s 

variation  in  shock  position  may  be  of  order  unity;  the  shock  can  move 

completely  through  the  channel,  as  opposed  to  case  (1)  where  the  shock 

motion  is  a perturbation  about  a steady  state  location,  x^.  Here,  for 

convenience,  we  again  write  the  steady  state  shock  wave  location  as  x^ 

and  write  for  the  general  shock  location, 

X = x +e'x,+...  =x„+A.x  +Fx,  + ...  (15) 

s so  si  U so  si 

Again,  (he  y dependence  of  the  shock  shape  is  of  higher  orrier  than  that 

retained  in  cqn.  (15). 

Because'  tlie  outer  solutions  do  not  satisfy  the  shock  jump  condi- 
tions to  second  order,  it  is  necessary  to  considi-r  an  inni-r  region  en- 
closing the  shock  wave,  just  as  in  case  (1  ).  ^ However,  ()io  situation 
in  this  case  is  complicated  by  the  fact  that  the  shock  wave  moves  over 


distances  of  order  unity.  That  is,  in  case  T),  liecause  the  shock  motion 


is  ,»  perturbation  alioul  \^y  witli  the  f)r<li'r  of  the  distance  moved  small 
compared  to  tt\e  lliickness  of  tlie  inner  re^>ion,  it  is  possible  to  con- 
sider a stationary  inner  remon,  inside  which  lliere  is  a higher  order 
stiock  w.ive  motion.  Here,  on  llie  other  hand,  one  must  consider  an 
inner  region  (ui  a coordinate  system  which  movies  at  the  shock  wave 
velocity;  in  view  of  the  I'xpansion  for  x^,  eijn  (1'’),  it  is  sufficient  to 
considei-  a .oordinale  system  moving  at  lh>-  vi  locity  dx^^/d  1 . 

If  (piantities  relative  to  tlie  moving  coordinate  system  are  de  - 

1 t t '1  ba , b , c , d I 

^ V ! 1 hi- . 1 , g ) 


noted  by  a cai'id,  then, 
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<4  - 1 X 

’ so 
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h . 


h - u X 
\ so 


so 
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where  p is  the  vector  vlocily  and  h^  the  st..gn.ition  enthalpy  in  the  mov- 
ing coordinate  system;  h^  h ' c/'/i,  and  q ..  re  the  stagnation  enthalpy 
ami  M'locitv  in  the  absolute  systems  respectively,  ami  h is  the  static 
enthalpy. 

Since  thi'  moving  system  is  a linearly  ,u  celebrating  system,  the 
conservation  eijuations  may  lie  derived  from  tliose  valid  tn  .in  im'ftial 
system  by  using  the  simple  transformations. 


1.  < --  --  — — 17a,  b,c) 

i)  1 P I'  so  ()  s;  dx  dy  Py 

Ihus,  the  mass,  momentum,  and  ene  r gy  c on  s e r valion  equations  becomi 


I) 


llH.i  I 


3T 


(q  • V)q 


1 '' 

- — V - i X 
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dh 
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A 

9T 


+ q 


''-A  .. 

^ h + u X 
t so 
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YP  9T 


(18b) 


(18c) 


Also,  the  following  dimensionless  property  relations  hold  for  a perfect 
gas  with  constant  specific  heats: 

T 


/-I 


]_ 

~ Y*1 

p = T ' e 


(19a  I 


(19b) 


P = P T (19c) 

where  s is  the  specific  entropy  made  dimensionless  with  respect  to  R, 
the  gas  constant.  If  we  now  write  the  velocity  as 

(20) 


-►  A A ^ 

q - 


Then,  using  eqns.  (18c)  and  (19),  one  can  show  that  the  energy  e quation 
can  be  written  as 


q • V(h,  + 'M  + X X ) 
^ t T so 


P 9 As  ^ a 


YP  9T 


- q 


9T 


(21) 


where  ux  ^-ix  -q-Vxx  . Now,  the  Reynolds  numl>er  is 
■ so  so  so 

tal<en  to  be  large  enough  that  viscous  effects  are  negligible,  and  the 

flow  is  transonic  so  that  shock  waves  are  weak  (i.e.  , As  = 0(r  )).  If 

one  writes  the  equation  expressing  the  change  in  vorticity  along  a 

streamline  in  the  moving  coordinate  system,  one  can  thus  show  ill  at 
7 ^ 

q = 0(r  ),  Moreover,  since  t “ 0(c  ),  it  is  easy  to  show  that  the 

a ^ 

terms  on  the  right  hand  side  of  eqn.  (21)  are  negligible;  that  is  to  the 


IS 


order  desired,  eqn.  (Zl)  hci  omes. 


h f 4) A I X X = F(T)  (id 

I I so 

Ti'is  is  the  form  of  the  Bernoulli  equation  to  be  used  in  the  moving 
coordinate  system.  The  corresponding  equation  in  the  inertial  frame 


hj  + = F(T)  Ui) 

A A 

The  relationsliip  between  F(T)  and  F(T)  is  easily  found  by  writing 
equation  {id)  in  terms  of  absolute  quantities,  using  equations  (16),  (17t, 
and  the  relation  between  4>  and  4>,  i.e. 

. A 

4>  = 4>  + X .X  (d4  ( 

s o 

and  by  replacing  h^^  +■  by  F(T)  according  to  eqn.  (23).  The  result  is. 


F(T)  ♦ 


= F(T) 


.2  2 4 ,2  4 A 

Actually,  since  x = k c (d  x /dt)  = 0(r  ),  F'(T)  F(1  ) to  the  order 

s o 2 so 

f A 

desired.  The  values  of  l''(T)  and  F(T)  upstream  of  the  shock  wave  can 
be  found  easily  by  noting  that  there  thj  flow  is  steady,  so  0,  and 

F = constant  evaluated  in  the  undisturbed  flow.  Thus 


F (T)  F (T)  - (261 

u u 2 2(y-l) 

A 

The  values  of  F and  F downstream  of  the  wave  may  be  determined  by 
noting  that  the  stagnation  enthalpy  relative  to  tlie  shock  wave,  written 


^Guderlcy,  K.  G.  , The  Theory  of  Transonic  Flow,  Pergamon  I'ress, 
Addison- Wesley,  1962,  p.  7. 


lUtlier  in  innex'  or  outer  variables,  is  conserved  across  the  wave,  i.e. 

■\  A 

h , = h . where  the  subscripts  u and  d refer  to  positions  immediately 
upstream  and  immediately  downstream  of  the  wave,  respectively. 
Then,  using  eqns.  (16),  one  can  show  that,  therefore. 


X. 


. I 

X 

so  , • , s 

h,  , - U X + — r h - U X + — 

Id  d so  c.  tu  u so  2 


(27) 


and  if,  again,  we  let  C = u,  (x  ) andnote  that  from  the  first  order 

u 1 u so 

solutions,  across  the  shock,  u = -u  , then  eqn.  (27)  becomes. 

Id  1 u 


> - 2e  X C t-  . . . 
td  2(y-l)  so  u 


(28) 


where  h = — ~ and  u - I + e u,  + . . . have  been  used.  Now,  since 

t u 2 ( Y - I ) 1 

the  unsteady  term.s  in  the  outer,  inertial  frame,  solutions  are  of  sec- 

2 

ond  order,  and  = 0(r  ),  it  is  seen  that  downstream  of  the  wave, 

4 3.2, 

4'._  = 0(e  ).  Hence,  to  order  c (since  x = k_,<  d x /dt),  one  finds 
T so  2 so 


from  eqns.  (23)  and  (25)  that 


. 2 

A f c O 

F.(T)  --  F.(T)  - -p 
d d 2 2 (y- 1 ) 


-liL 


- 2e  X C + . . . 
so  u 


(29) 


H”nce,  the  Bernoulli  equations  in  the  moving  coordinate  system,  eqn. 
(22),  and  in  the  stationary  coordinate  system,  eqn.  (2  1),  are  known, 
Ijoth  upstreatn  and  downstream  of  the  shock  wav(?,  by  virtue  of  eqns. 
(26)  and  (29).  It  should  be  noted  that  the  same  results  are  found  if 
detailed  inner  region  calculations  arc  carried  out. 

An  inner  region,  enclosing  the  shock  wave,  must  Ije  etjnsidered 


because,  as  mentioned  previc)usly,  the  outer  sedutions  do  not  satisfy 
the  shock  jump  conditions  in  second  order.  In  this  thin  inner  region, 


then,  the  solutions  must  (i)  match  with  the  outer  solutions  term  by- 
term  as  the  outer  region  is  approached  in  an  appropriate  limit,  tii) 
satisfy  the  wall  boundary  conditions,  and  (iii)  satisfy  the  shock  jump 
condition  at  the  shock  wave.  The  analysis  differs  from  that  followed 
in  reference  (1)  only  in  a few  minor  instances,  so  a very  brief  review 
is  given  here,  with  emphasis  on  the  those  differences. 

The  gas  dynamic  equation  in  the  moving  coordinate  system  may 
be  derived  from  eqns.  (18),  (19),  and  (20),  with  q =0  and  As  -■  0 

3. 

(to  the  order  desired  here).  However,  it  is  simpler  to  transform  the 

(4  ) 

equation  written  in  terrns  of  absolute  quantities,  using  eqns.  (17)  and 
(24).  The  result  is, 

Z ^ A^A  A 

(a  - + (a  - - 2$.  <*>.  ■ 

X XX  y yy  ^ y r i 

, (30) 

- 24>/,$a*  - 2$A$n4,  - X - X X + X X =0 

X xT  y yT  so  x so  so  so 

Just  as  in  case  (1)^^^,  it  can  be  shown  that  the  inner  region  is  of 
order  ( ^ in  thickness.  Hence,  inner  variables  are  written  as  follows: 


(x-x  )c 
s o 


= ^2'  ■^2' 


(1  la,  b, c) 


and  an  inner  velocity  potential  may  be  defined  as 

A*:  - 1 / 2 ^ A - . -1/2 

<I'(x,y,T)=(  <I>(x,y,T)-x(l-x  )+e  x 

so  so 

;;;  ^ ^ ^ 
-t-  4>  (x  , y , t ) 

whi?re  perturbation  potential.  From  eqns.  (.5  1)  and  (3Z),  it 

is  seen  that 


u = '1>.\  - 4-  = u = 1 t (j) 


* so 


( i ia) 


^ l/Z*  1/2  =:=  1/2 

4>,'  - t <t>  r:  r V - r , 

y yV  y 


( 1 •>  b) 


The  matching  conditions  for  u'  and  v , i.e.  the  solutions  to 
which  they  should  match  as  x'  — oo,  are  found  by  transforming  the 
outer  solutions  for  u and  v (eqns.  (4),  (5a),  (5b),  and  (5c)  with  ^ - 0) 

into  velocities  in  the  moving  coordinate  system,  u and  v 'eqns.  flffl 
and  (I6g)),  then  expanding  these  solutions  about  x^^  and  writing  the 
resulting  expressions  in  inner  variables.  If  x^^  is  subtracted  from 
eqn.  (ilia)  in  reference  (1),  then  eqns.  (31)  in  that  reference  are  pre- 
cisely the  desired  relations.  Now  in  case  (1),  because  x^  - constant 


th 


w 


as  the  zero  order  shock  wave  position,  no  further  matching  was 


required.  Here,  because  x = x (t),  there  is  a function  of  time  in 

^ SOSO 

the  potential  which  must  be  found,  in  order  to  evaluate  properly 

in  the  inner  region  counterpart  of  eqn.  (22),  the  Bernoulli  equation. 

This  term  is  found  by  matching  the  potential;  i.e.  using  eqns.  (24) 

2 

and(32),  with  4>  - x + ^ 4>2  • • • » expanding  "t*  about  x^^,  one 

* 

can  show  that  the  inner  perturbation  potential,  4>  > niust  approach 


>:-■  * « 1/2 
(1  oo/  y . f ) = c (x^^)  )-  . . . 


(34) 


where  6i  value  downstream  of  the  shock  for  the  upper  sign 

1 so 

in  (()  , and  its  value  upstream  of  the  shock  for  the  low’  r sign.  F'inally, 


the  resulting  expansion  for  c()  is, 

* 1/2  * 3/2  * , 2 

4.  r +e4>i  ^ ‘■‘’2  ^ 


(351 


where,  since  x = x (t),  then  4k  (x  ) is  a function  of  time  alone, 
so  so  Iso 


and. 


^ * so 

•P  * 
t ' 


4’,  ). 

d I lx  so 
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with  d)  (x  ) = C upstream  of  the  shock  and  4>  (x.  ) = - C down- 

Ixso  u Ixso  u 

stream  of  the  shock.  It  should  be  noted  that  in  case  (2.),  C = C (t) 

u u 

since  .x  = x (t). 
so  so 

5jc 

With  the  e.xpansion  for  known,  the  governing  equation  in  the 
inner  region  may  he  derived.  Thus,  eqn.  (22),  with 


h = h + 
t 


,-2  o2  2 a2  /-2 

u +v  _ a ^ u +v 


is  used  for  a in  eqn.  (30),  and  the  resulting  equation,  using  eqns.  ib) 
(31),  (32),  and  (35),  is  written  in  terms  of  inner  variables;  finally, 
governing  equations  for  each  4*^^  may  be  written.  These  equations 
are  precisely  the  same  as  those  derived  for  case  (1),  because  ev'en 
though  u,  (x  ) - C is  not  a constant,  and  in  fact  is  a function  of  time, 
there  arc  no  time  derivatives  involved  in  the  governing  equations  to 
the  order  desired,  as  a result  of  the  order  of  t^.  Also,  the  boundary 
conditions  at  the  wall  and  the  shock  jump  conditions  are  unchanged. 
Hence,  the  inner  region  solutions  are  the  same  as  those  given  in  ref- 
erence (1);  that  is,  one  can  define  an  added  second  order  potential 
function,  , in  the  same  way  and  obtain  for  it  the  same  solution. 
Therefore,  eqns.  (4),  and  (5a)  - (5f)  give  a composite  general  solution 
for  u and  v,  uniformly  valid  to  second  order,  for  case  (2)  as  well  as 
case  (1).  There  remains  only  the  question  of  the  proper  forms  for  h 
and  the  shock  position,  corresponding  to  eqns.  (5f)  and  (7a)  respect- 
ively, in  case  (2 ). 

The  equation  for  h , which  occurs  in  the  seiond  order  outer 

X 

solution  for  u,  is  found  by  deriving  the  third  oi'dcr  velocity  solutions 


and  applying  the  boundary  conditions.  The  procedure  follows  exai  tly 
those  steps  given  in  reference  (1)  and  so  will  not  be  repeated  here. 

The  equation  found  for  h^  is,  then,  for  case  (Z), 

h = -■-({"  I-  (Zy-  3)u^)  ^ -AILL  (37) 

X 6 ’ 1 u^ 

where  A(t)  is  a func  tion  of  integration.  This  is  the  same  equation  found 
for  h in  case  except  that  in  case  (1),  an  additional  term,  h , ap- 


pears  in  the  equation.  Because  1 


T 1'  Dh/OT  — (9h/t)tl,  it  is 


I Tear  that  in  case  (Z)  the  term  9h/8t  appears  in  a higher  order  calcula- 
tion rather  than  in  the  same  equation  as  h^;  hence  in  case  (2),  the  cal- 
culation for  h is  easier  than  in  case  (1).  The  boundary  condition  to  be 

X 

applied,  in  order  to  evaluate  A(t),  is  that  a given,  unsteady,  pressure 
is  impressed  downstream  of  the  flow,  say  at  x - X,  these  pressure  os- 
cillations being  in  second  order,  as  stated  earlier.  Now,  since  one 
can  write  the  pressure  as 

and  since  the  unsteady  part  of  u^  is  lt^(x,  t),  it  is  clear  that  a condition 
on  h is  equivalent  to  a condition  on  the  second  order  pressure;  hence 

X 

boundary  conditions  on  h^  are  given  downstream  of  the  wave. 

For  steady  flow,  A(t)  is  a constant,  say  c^,  so 

C 

h = - -^(f"  ^ (2y-3)uf)  t — (39) 

X b 1 Uj 

Thus,  eqn.  (39)  holds  upstream  where  the  flow  is  steady,  and  whore 

C ■ C , sav.  Downstream  of  the  shock,  wliere  the  Ilow  may  be 
2 2u 


I 

J 


unsteady,  it  is  c:onvenient  to  write  as  the  sum  of  a steady  term, 

(h  I . eiven  by  eqn.  and  an  unsteady  contribution,  li  , whc?re, 

X s s ^ 


A(t)-C 


h (X  , t)  - 
X 


If  the  boundary  condition  is  written  as 


h (X,  t) 

X 


Uj(X) 


1 hen  it  is  seen  that  Ad)  - - G(t)  and  so  in  general, 

h (2y-3)u  J t-  — - t-  (4£) 

X 6 1 u j u j 

When  this  solution  is  compared  with  the  corresponding  solution  for  case 
(1),  given  t)y  eqn.  (5f),  it  is  seen  that  the  only  difference  is  the-  absence 
of  i)i«^  time  lag  in  G,  in  eqn.  (4<d).  This,  of  course,  arises  as  a result 


of  the  fact  that  in  the  present  case  T ^ that  in  the  limit  as 

r -*  0,  there  is  no  time  lag;  the  impressed  pressure  oscillations  and 
the  corresponding  oscillations  in  u,  f’,  p , etc.,  are  in  phase. 

The  equation  for  the  shock  wave  location  is  found  by  applying 
the  mass  conservation  principle  to  a control  volume  cont. lining  the  shock 
wave,  and  hence  attached  to  the  moving  <-oordinate  system.  If  this  con- 
trol volume  is  considered  to  liave  an  extent,  in  the  flow  direction,  of 


order  f (th(!  control  volume  e.xtends  from  a constant  negativ'c  valu(' 

of  ,x  to  a constant  p(i)itive  value  of  x ),  then  the  variation  with  time  of 


the  ma.ss  within  the'  control  volume  is  at  least  of  oidei  f 


.Since  the 


interest  here  is  only  in  terms  to  and  including  0(r  ),  the  time  vaiiation 


IL.- 


k 


u-rm  need  not  be  considered,  and  so  conservation  of  mass  demands 
that 


\v  , A 

l>  It  d y 


'V  * , A 
f>  u d y 


lA  i I 


where,  from  eqns.  (4),  (5),  (16),  (19),  (d^),  (^9),  (3d),  and  (i(.). 


p u = 1 - f 


M,  r 

u SO  6 1 


t (YH)Uj^f"^  ^ \ =:=  • '^z'^soy  ^ 


V 

I dsoj  • • • 

Upstream  of  the  wave,  u,  = C and  C ^ = 0,  while  downstream  of  the 

In  * 


.)■ 


(44  ) 


wave,  u - - C and  ? ^ is  calculated  using  eqn.  (Sd)  and  evaluated  at 


some  X - constant.  In  both  locations,  f"  = f"(x  ) + . . . = f"  , say 

so  ^ ' 


so 


Finally  As  = s - s =0  upstream  of  the  wave  and  is  equal  to  the  As 


(5) 


across  a weak  shock  downstream  of  the  wave,  i.  e. 


,-s  ) = c — ^ (y+1  )C^  + . . . 
d u i ' u 


(45) 


Since  we  are  interested  only  in  terms  up  to  O(r^),  and  y r | + O(c^), 

w 


it  is  seen  from  eqns.  (4  3)  and  (44)  that  in  the  integral  in  eqn.  (4  1), 
the  limit  y^  may  be  replaced  by  1.  Now,  if  eqn.  (44)  is  substituU‘d 


into  eqn.  (4  3).  using  the  above  mentioned  sul)stilutions  upstream  and 


downstream  of  the  wave,  then  from  the  first  nonzt^ro  term  (0(f  )),  it 
is  found  that  since  f ^ ,,.dy  - 0, 

t)  X 


Fiepmann,  H.  W . and  Roshko,  A.,  F:ietnents  of  Ciasdynamics , 
John  Wiley  and  Sons,  1957,  p.  60. 


Z3 


I 

J 


4k. 


dx 


(Y+1)  dt 


= h + h 

X X 

d u 


.5  U 


(46) 


Thus,  the  form  is  exactly  that  found  for  case  (1),  eqn.  (7a).  However, 
in  eqn.  (46).  C and  f"  , for  example,  are  functions  of  time  rather 
than  being  constants  as  in  eqn.  (7a);  as  will  be  seen,  this  leads  to  dif- 
ferences in  shock  wave  velocity  between  the  two  cases. 

It  should  be  noted  at  this  point  that  the  general  form  of  the  solu- 
tion is  the  same  for  the  two  cases.  'I  he  differences,  which  are  funda- 
mental, are  in  the  details  of  the  solution,  i.e.,  in  case  1)  a time  lag, 
t^  , is  involved  and  in  case  (2)  it  is  not,  and  in  case  (1)  functions  evalu- 
ated at  the  shock  are  constant,  while  in  case  (2)  they  are  functions  of 
time . 


Before  showing  example  numerical  calculations  for  case  (2),  it 
is  of  interest  to  write  equation  (4  6)  in  a more  convenient  form  and  com- 
pare it  with  its  counterpart  from  case  (1).  The  problem  c onsidered 
is,  again,  that  of  a flow  which  accelerates  through  sonic  velocity  at  a 
throat.  There  is  a shock  wave  in  the  flow,  which  is  steady  until,  at 
time  t = 0,  second  order  pressure  oscillations  are  impressed  upon  the 
flow  downstream  of  the  shock  wave.  If  eqn.  (42)  is  used  to  calculate 

h and  h , then,  since  C,  =0,  eqn.  (46)  becomes, 
d u 


4k,,  dx 
2 so 


(ytl)  dt 


C,  f + G(t) 

2 , 5 u 


(4  7) 


Equations  (5f)  and  (7a)  may  be  used  to  derive  the  corresponding  equa- 
tion for  case  G ).  Thus, 


24 


4 k , dx  , 
1 s I 

( y+ 1 ) d t 


C Z , i uo  to 

uo  1_  d 


where  t,  = t.(x  ) (see  eqn.  (5g)),  and  where  C - C (x  ).  At  steady 
to  t o , , uouo 

dx  dx 

state,  when  G = 0 and  — ~ 0 = — both  eqns.  (4  7)  and  (48) 

dt  dt 


Sive, 


C,  = - C^x  ) = - ^ (4« 

Z , 3 u o 3 uo 

d 

so  that  finally,  one  can  write,  for  eqns.  (47)  and  (48),  where,  again, 

C = u,  (x  ), 
u Iso 


2\  ^3 


4k^  d 

(Y+*)  d 


so  _ 1_  I 

t C \ 

u ' 


4k,  dx 
1_  s 1 

(Y+1)  dt 


^ (C^  - - G(t) 

3 uo  u , 


uO 


(50a) 


(50b) 


The  differences  between  the  two  shock  motions  are  apparent,  when 
written  in  this  form.  A more  detailed  discussion  of  the  shock  motion 
for  (,'ase  (Z),  eqn.  (47)  or  eqn.  (50a),  is  given  later. 

Numerical  solutions  for  a typical  case  (Z)  flow  are  shown  in  fig- 
ures )a  (t'O,  steady  state),  6a  to  6c,  and  7.  The  wall  shape,  functional 
form  of  G(t),and  in  fact  all  parameters  except  for  t arc  the  same  as 
those  used  in  the  niimerical  examples  shown  for  case  |1),  fig\ires  3 to 
5 (Eqns.  9).  Thus,  in  Case  (Z),  = 100  while  in  case  (1), 

Tj  “ ZO.  The  most  significant  difference  noted  l)etwtM-n  (tie  (wo  castus 


is  the  overall  shock  wav(t  motion. 


Unified  Solutions  - Case  (i) 


[ 


Tile  differences  bi-tween  the  two  cases  considered  are  very  clear 
from  an  asymptotic  viewpoint.  Hcwever,  if  one  wishes  to  make  flow 
field  computations  for  a given  set  of  physical  constants,  the  choice  of 
which  set  of  solutions  to  use  is  not  so  clear.  Thus,  if  r 0.  1 and 
T - 40,  then  should  one  use  case  '1)  solutions  with  3 ^ . 2B  or  case 

(<d)  solutions  with  5?  A relatively  simple  answer  to  this  ques- 

tion consists  of  using  the  results  of  the  previous  analysis  to  write  a 
unified  solution  for  numcrii  al  computations. 

From  the  preceding  analysis,  it  is  clear  that  the  general  form  of 
the  solution  for  case  (1)  i.i  the  same  as  that  for  case  (^).  In  what  fol- 
lows, the  aim  is  to  write  a solution  which  reduces  properly  to  either 
case.  We  first  rewrite  the  time  lag,  eqn.  (^g),  as  follows: 


f —L.  dg 
c J (y+1) 

X 


(5-1  1 


Then  as  t -►  t ^ = (k^r)  , eqn.  (51)  reduces  to  eqn.  (5g),  and  as 

d -1 

r -►  T (k.f  ) , t 0(r),  and  becomes  negligibly  small.  '1  he  uni- 

c c.  f 

fied  solution  thus  is  writt»*n  with  Cift-t.)  in  h , wdth  t.  giv<'n  bv  eqn. 

P X P 


(51). 

The  only  remaining  differenci*  is  that  in  the  t:quations  for  1h<' 
shockwave  v<*locily,  i.e.,  eqns.  (4  7)  and  (48).  For  the  unified  solu- 
tion, we  write 


X 

s 


X,,  '<  X (t) 

U S 


(h^a) 





1 


4 


1 


dx 


r c 


2 (y+1)  dt 


— 

i C 

’ U 


2d  j u f o 


f o 


l_  __J_  _di_ 

T c J (yi-l)  u ) 


(52c) 


and  for  the  specific  problem  considered  here,  where  oscillations  are 

impressed  on  a flow  with  a steady  state  shock  location, 

C-,  = - 2vC^  / i and  so  eqn.  (52b)  becomes : 

2d  ‘ uo 


. dx 
4 s 


T e (y+1)  dt 


- C^)  - G(t  - t 


UO 


io ' 


(53) 


Thus,  for  T-^  T ^ - (kj  r)  , the  solution  of  eqn.  (53)  or  (5Zb)  is  for 

X /c,  i.  e.  , X ■-  c X , . Then  C = C + 0(e),  and  eqn.  (53 ) reduces  to 
“si 


uo 


eqn.  (SOlj).  On  tlie  other  hand,  if  t -*  - k^r  , t^  = 0(e)  and  eqn.  (53) 

t 

reduces  to  eqn.  (50a),  with  x = 0(1).  Therefore-,  the  unified  solution 

s 

sufiK^'-'^h^d  he  ri!  is  f^iven  byeqns.  (4),  (5a)  - (5f),  (51),  and  (52),  witli 
eqn.  (53)  replacing  eqn.  (52b)  for  the  specific  problem  considi-red  here. 

An  indication  of  the  relationship  between  calculations  carried  out 
using  case  i 1 ),  case  (2),  and  case  (3)  solutions,  is  given  in  figures  8 
and  9.  In  figure  8,  case  (l)(eqns.  (6)  and  (50b))  and  i-  ase  ( 3 ) ( eqns . 52  ) 
computations  for  shock  position  are  compared  for  t 20,  r - 0.0“ 

(i.  e.  , kj  1 in  case  (1)).  In  figure  9,  case  (2)  and  case  (5)  computa- 
tions for  shock  position  and  shock  velocity  are  compari'd  ff)  r r 400, 
c 0.05  (i.e.,  k^  - 1 in  case  (2)).  In  figuia-  9,  the  difference  in  shock 

velocities  is  entirely  due  to  the  difference  lu-twei'n  G(t)  anrl  G(t-t,  ) 

' to 


and  thus  gives  an  indication  of  the  effect  of  carrying  t,  in  G in  case  (3), 

f 


;» 

I 

1 


1 


i_e.  , t ^ t In  cai'li  l om pa ri.son,  case  (3^  is  considered  to  be  the 

to  ~ ( 

more  accurate  calculation,  which  reduces  to  the  case  with  which  it  is 
being  compared,  as  r -*  0.  It  should  be  noted  that  as  t increases  in 
value  to  the  point  that  the  difference  in  u,  say,  between  values  found  using 
in  G and  not  using  in  G,  becomes  of  order  e , then  t^  should  be 
neglected;  i.e,  in  the  solutions,  terms  of  order  have  been  neglec- 
ted, so  a term  of  this  order  should  not  be  carried  as  a correction  to 
G(t). 

In  summary,  then,  it  is  possible  to  construct  a unified  computa- 
tional sche^me  from  which  numerical  results  may  t)c  oiilai;ied  for  a 
large  range;  of  ( and  t values.  Thi-se  solutions  are  valid  in  the  so 
called  slowly  varying  time  regime,  which  covers  a range  of  character- 
istic limes,  for  impressed  pressure  oscillations,  of  cons  ido  raljle 
technic^ll  interest. 


Large  Amplitude  Shock  Wave  Motion 

As  indicated  previously,  when  t = 0(f,  ) (case  the  shock 

motion  resulting  from  pressure  oscillations  impresst'd  downstream  of 
the  shock  wave,  has  an  amplitude  of  order  unity.  As  a r.  ■ ull,  there 
are  coiiditions  under  wtuch  th«'  shock  will  move  upstrc'atii  through  the 
noz/.le,  disappear  upstream,  and  tht'n  reappear  as  the  flownstream 
plenum  pressure  drops  to  the  point  where  a shock  wavt-  in  the  channel 
is  necessary  to  satisfy  this  instantaneous  pressure  requi  ■ ••ment . ’1  he 

i-onflitions  for  this  oceurrenee  and  the  subs<'quent  shock  wavi-  motion 
depenrl  in  a complex  manner  upon  (he  amplitud(’  of  (he  for.  ing  hinctiiui, 


ZH 


G,  the  steady  state  condition  about  which  the  oscillations  occur,  repre- 
sented by  the  wall  shape,  f(x),  and  the  numerical  value  of  the 

time  constant,  represented  by 

The  equation  which  governs  the  shock  motion,  repeated  here  for 
convenience,  is  eqn.  (47),  or  in  the  form  suited  to  the  proljlem  under 
consideration,  eqn.  (50a). 


' fc  - ' call 

rvTn  ' C L 2d  i >.  _ 


4 k dx 

^ 

(y^I)  dt 


— = ^ (C^  -C^  - G(t) 

C 3 uo  u 


(50a) 


where,  since  .x  = x + 0(e),  then  to  the  order  considered  here,  x and 
s so  s 

X are  interchangeable.  The  interesting  point  indicated  by  the  e.xis- 
so 

tence  of  these  equations  is  that  in  spite  of  the  fact  that  signals  from  the 
impressed  disturbances  reach  the  shock  wave  "instantaneously,  " that 
is  there  is  no  lag  time  in  the  pressure  or  velocity  solutions  for  e.xarn- 
ple,  the  shock  wave  does  not  respond  instantaneously.  The  shock  ve- 
locity is  finite,  and  indeed  there  is  a lag  between  the  impressed  dis- 
turbance, G(t),  and  the  resulting  stretched  shock  velocity  dx  /dt,  due 

s o 

to  the  other  terms  in  the  equation.  Thus,  in  equation  (50a),  for  exam- 
ple, it  is  seen  that  the  term  (C^  - C^)  always  has  a sign  such  that 

" uo  u 

its  effect  is  to  cause  the  shock  to  move  toward  the  equilibrium  or 
steady  state  position.  On  the  other  hand,  G(t)  is  a forcing  function 
which  changes  sign  periodically.  The  result  is  a shock  motion  which 


lag.s  G(t). 


It  is  clear  from  eqn.  (47)  or  (SOa)  that  sinj>ularities  occur  as  the 


shoe  k wave  approaches  the  throat,  and  - 0.  The  behavior  of  inte- 
gral e urves  which  cross  the  axis,  x = 0,  on  the  x - t plane  can 

so  s o 

be  found  for  a sinusoidal  Gft)  by  writing  eqn.  (4  7)  for  x (and  thus  C ) 

so  u 

small  compared  to  unity  and  for  t - t « 1.  where  t is  the  value  of  t at 

o o 

which  the  bracket  on  the  right  hand  side  of  eqn.  (47)  goes  to  zero  at  the 
throat,  0.  'Ihus,  if,  for  example. 


then 


G G sin  bl 
o 


sin  bt  :r  - C , ,/G 
o 2d  o 


and  eqn.  (47  ) becomes,  for  x « 1 and  t - t « 1 

so  o’ 


dx 


so 


dt 


(\  '■1)  ,, 

-r. ( G c o s b t 

4 k^  o 


(hi) 


( 55 1 


(56a) 


. o»  1.  |56bl 

where,  again,  C u (x  t is  the  value  of  u at  x , upstream  of  the 

II  lUSO  1 

shock,  and  where  « t - t^.  In  the  neighborhood  of  the  throat,  a 
parabolic  wall  shape  is  representative  of  actual  practice  Thus,  a 
typical  wall  shape  and  the  corresponding  solution  to  eqn.  (56)  are,  in 
the  neighl)orhood  of  the  throat. 


Thus  if  cos  hi  > t),thi-  point  0,1  is  a centcT  and  the  integral  curves 
o o 

(ellipses)  in  the  neighborhood  of  this  center  cross  ~ 0 with  an  in- 

finite slope.  Such  curves  are  shown  in  figure  10.  On  the  other  hand, 
if  cos  bt  < 0,  (he  integral  curves  in  the  neighborhood  of  0,  t^  are 
hyperbolae,  the  point  being  a saddle  point,  as  illustrated  in  figure  10, 


and  the  two  integral  curves  pass  through  the  point  0,  t^  with  slopes 


ft  - -\^5/2 


3/Z  b G ) 1/Z 


/2,  1/2 
a 


An  understanding  of  the  possible  shock  motions  may  be  gained  by 

analyzing  the  integral  curves  which  pass  through  the  saddle  points.  The 

three  possible  configurations  for  these  curves  are  sketched  in  figure  11. 

In  these  sketches,  the  arrows  indicate  the  direction  the  solutions  mvist 

follow  as  time  increases.  In  figure  11a,  conditions  are  such  that  the 

integral  curves  entering  the  saddle  point  originate  from  a particular 

X = X at  t = 0.  Those  leaving  the  saddle  begin  to  rise,  then  reverse 
so  o 

their  directions  and  cross  the  time  axis  with  vortical  slope  at  some 
point  between  the  center  and  Ihe  next  saddle  point.  Other  integral 
I'urves  are  sketched  also,  in  dotted  lines.  As  indicated  in  the  sketch, 
the  paths  traced  by  the  integral  curves  are  repetitive.  In  figure  Me, 
the  opposite  situation  exists;  the  integral  curves  entering  tlie  saddle 
point  begin  on  tlie  abscissa,  lietween  a saddle  and  a center,  with  an 


infinite  slope  and  then  change  direction  and  enter  the  next  saddle  point. 
Those  curves  li-aving  the  saddle?  never  return  to  the  axis  x^^  0,  Init 

asymptotii-ally  approach  a single  periodic  curve  (for  a given  C, 


1 


Ihosi'  i-urves  which  orifiinate  with  an  x greater  than  any  x on  this 

so  so 

piTiocJic  curve  will  approach  the  periodic  curve  asymptotically  from 

above.  This  periodic  curve  is  nearly  symmetric  about  x^,  the  steady 

statt'  value  of  In  the  dividing  case,  shown  in  figure  lib,  the 

curves  entering  and  leaving  the  saddle  points  are  the  same  curve. 

The  integral  curve  jnap  obtained  in  any  given  case  depends  upon 

the  values  of  the  amplitude  of  G,  C , k , and  the  value  of  dC  /d.x 

lid  c u s o 

Although  general  solutions  from  which  a general  criterion  for  the  di- 
viding condition  (figure  lib)  could  be  derived  are  not  available,  an 
approximate  result  can  be  found  for  G as  given  in  eqn.  (54)  and  f(x) 
as  in  eqn.  (57a).  Thus,  with  these  substitutions,  eqn.  (4  7)  becomes. 


K X 


dx 
so  dt 


- - C - r X -G  sinbt 
cd  so  o 


(59a) 


K 


,5/^ 

2 


( V+ 1 ) 


5/2 


(59b, c) 


and  the  slopes  of  the  integral  curves  at  the  saddle  point  are  given  by 

eqn.  (58).  It  is  assumed  that  the  integral  curve  which  passes  through 

the  saddle  point  at  bt  - bt^  and  also  through  the  next  saddle  point  at 

^ bt^  +-  Z 7T  (e.g.,  see  figure  lib),  is  approximately  symmetric  about 

bt  ” bt  *■  7T.  At  this  point,  x is  taken  to  have  its  maximum  value, 

(x  ) . and  so  dx  /dt  = 0.  Hence,  from  eqn.  (59a), 
so  m so  M ' /> 


(x  ) 
so  m 


1/j 


; 


. 
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(60) 


The  integral  curve  in  question  is  approximated  by  a polynomial. 


1 


~ ~ ~2  ~ 2 
X - C.t(l-t)  I 4C,t  (l-t) 


f6  I a ) 


I b(t-t^)/2;r 


(bib) 


which  satisfies  the  conditions  that  x=0  at  bl  - bt^  and  bt  - bt^^  f 2 r,  and 
dx  /dl  -■  0 at  bt  = bt  + n.  Moreover,  at  bt  = bl  and  at  bl  bt  t In, 

^ r\  O O 


SO 


the  correct  slopes  of  the  integral  curves  are  obtained  if  bC^/27r  is 

given  by  the  right  hand  side  of  eqn.  (58)  with  the  upper  sign.  1-inally, 

at  bt  = bt  + .7, 
o 


X =(C,+C)/4  = (x  ) 
so  1 2 so  m 


Since  C,  and  (x  ) are  known  in  terms  of  the  desired  parameters, 
1 so  m 

it  is  nece.ssary  now  to  find  C^,  or  W'hat  is  simpler  C^/C^,  where 


4(x  ) 

so  m 


~ 1 


(62) 


"'I  ^1 

The  desired  relation  is  found  by  integrating  eqn.  (59a),  over  the  period 

bt  to  bt  + Ztt.  The  result  is,  using  eqn.  (61b), 
o o 


/‘ 


X dt 
0 so 


2d 

r 


(63) 


If  eqn.  (61a)  is  usc-d  in  evaluating  eqn.  (63),  and  is  replaced, 

using  eqn.  (62),  the  following  cubic  equation  and  solution  are  found 

for  (x  ) /C  - R,  say, 
so  m 1 


(4B-1) 


i 
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(64a) 
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^1 


(64b) 


Krot>i  ccjn.  (6Z),  aftc'r  substituting  for  and  finds 

ihu  i-qu£ition  relating  ^ special  dividing  case 

illustrated  in  figure  lib. 

Z Z , \ 4/3 

q“  - C,  , t-  0.  41652 ( ) (65) 

^ (Y+l)a  V Y / 

1 hen  for  G greater  than  the  special  value  given  in  eqn.  (65),  the  inte- 
o 

gral  curves  shown  in  figure  11a  result,  while  for  less  than  the  spe- 
cial value,  thosi'  in  figure  11c  are  found. 

Plxample  calculations  of  the  int(*gral  curves  through  the  saddle 
points,  with  the  sinusoidal  forcing  function  given  in  eqn.  (54),  with 
parabolic  walls  as  in  eqn.  (57a),  and  using  the  approxinaate  form  to 
calculate  the  special  value  of  G^,  eqn.  (65),  are  shown  in  figure  12. 

The  calculations  were  carried  out  by  numerically  integrating  eqn. 
(59a),  using  eqn.  (58)  to  find  an  initial  condition  near  x^^  --  0.  In  the 
calculations,  b 2,  • a = (Y't^l)/2  - 1.2,  x^  = 1.5,  - (2f(x)/(Y+l )) 

and  C - 2vC‘  /i.  In  figure  12,  the  letters  a,  b,  and  c,  refer  to 

2d  ' uo 

the  corresponding  cases  sitown  in  figure  11.  In  eai:h  case,  only  the 

curves  through  one  saddle  point  are  shown;  the  repetitive  nature  ol 

the  curves  at  each  saddle  point  is  not  shown,  for  clarity.  It  should  be 

noted  that  the  value  of  t in  figure  12,  referring  to  the  location  of  a 

o 

saddle  point,  is  different  for  each  case.  1 he  centers,  which  also  oc- 
cur at  different  values  of  t for  each  case,  arc  noted  in  figure  12.  With 


the  parametric  values  given  above,  it  was  found  from  eqn.  (65)  that 
the  special  value  of  which  gives  the  behavior  shown  in  figure  lib. 


(G  ) = 4.968098  (66) 

o sp 

The  curve  labeled  b in  figure  12  indicates  that  this  value  is  quite  accu- 
rate. The  integral  curve  leaves  x =0,  t-t  0 and  returns  to 

so  o 

X = 0 very  nearly  at  t-t  - tt',  the  approximations  employed  in  de- 
so  so 

riving  eqn.  (65)  appear  to  be  justified.  The  curves  labeled  a and  c 

in  figure  12  were  calculated  using  G 5.  5 > (G  ) and  G 4 < (G  ) , 

o o sp  o o sp 

respectively.  In  each  of  these  cases,  curves  entering  and  leaving  the 
saddle  point  at  t-t  --  0 are  shown,  the  behavior  in  each  case  following 
that  sketched  in  the  corresponding  part  of  figure  11.  It  should  be  noted 
that  the  solutions  shown  in  figure  12  are  for  very  simple  (parabolic)  wall 
shapes,  and  that  no  simple  way  of  predicting  (61^)  for  the  case  of 
general  wall  shapes  exists.  However,  these  simple  calculations,  the 
results  of  which  are  shown  in  figure  12,  serve  to  illustrate  thi'  various 
cases  which  may  occur,  as  sketched  in  figure  11,  and  thus  are  extreme- 
ly useful  in  interpreting  results  found  for  more  complicated  geonnUries. 
For  such  geometries,  it  is  necessary  to  integrate  numerically  along  an 
integral  curve  leaving  a saddle  point  to  see  w'hich  case  occurs  for  ihi> 
given  parameters.  Examples  of  such  calculations  arc  shown  later. 

With  the  n^athematical  behavior  of  the  integral  curves  through 
saddle  points  understood,  it  is  possible  to  interpret  the  physical  behavior 


of  the  shock  wave  in  each  case.  Referring  to  figure  11a,  for  any  initial 
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I 
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condition  which  does  not  lie  on  an  integral  curve  entering  a saddle 
point  (two  are  illustrated  by  circles  in  figure  Ha),  the  shock  passes 
through  the  throat  and  disappears  upstream.  This  is  seen  by  following 
the  integral  curve  in  question  as  time  increases;  goes  to  zero,  for 

any  initial  condition,  between  a center  and  a saddle  point.  As  time 
increases,  then,  a saddle  point  occurs  at  - 0,  and  an  integral 
curve  rises  from  the  saddle  point  in  the  direction  of  t increasing. 

This  means  that  the  back  pressure  has  decreased  to  the  point  where 
a shock  wave  must  form  in  the  channel  in  order  to  satisfy  the  instan- 
taneous pressure  requirements.  This  is  seen  by  writing  the  pressure 
at  X = X,  using  eqns.  (38),  (5b),  and  (39);  since  f"(X)  = 0, 

2 |/’2v-3\  2 ^2d 

Pb  = ^ -^YUi(X)+e  y|(-^)uj(X)  - — 

where  Uj(X)  < 0.  From  eqn.  (67),  it  is  seen  that  the  conditions  for 
the  back  pressure  to  be  that  which  gives  the  subsonic  solution  (no  shock 
waves)  with  sonic  pressure  at  the  throat  (upper  dotted  curve  in  figure 
2 ) is 


l(X)  J 


c,  , + G(t)  = 0 


(68) 


But  this  condition,  for  the  case  where  G(t)  (and  hence  the  pressure)  is 
decreasing,  is  precisely  the  condition  for  the  saddle  point,  as  exempli- 
fied by  eqns.  (54)  and  (55)  and  the  discussion  following  these  equations. 
The  fact  that  this  back  pressure  requirement  must  be  satisfied  instan- 
taneously by  a shock  forming  at  the  throat  is  a result  of  the  fact  that 
there  is  no  time  lag  in  the  solutions  for  the  velocity,  pressure,  etc. 
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As  a result,  then,  the  proper  behavior  for  the  shock,  after  it  disap- 
pears, is  to  reappear  at  the  time  associated  with  the  first  saddle  point 
after  its  disappearance.  It  then  follows  the  path  jjiven  by  the  intc'gral 
curve  through  the  saddle  point  and  so  disappears  again,  forms  again 
at  the  throat  at  the  following  saddle,  etc.  Thus,  no  matter  what  thi- 
initial  condition  is,  the  rt'sulting  shock  motion  is  as.so<  iat<;d  witli  the 
integral  curves  leaving  the  saddle  points,  as  shown  in  figure  IJa.  lor 
the  periods  of  time  between  the  shocks  disappearing  upstream  and  re- 
forming at  the  throat,  the  flow  is  subsonic  throughout  the  channel.  If 
the  initial  condition  should  lie  on  an  integral  curve  entering  the  saddle 
point,  the  shock  moves  to  the  throat  and  moves  away  again  on  the  inte- 
gral curve  leaving  the  saddle  point.  Thereafter,  its  motion  is  tin.’ 
same  as  that  shown  in  figure  15a. 

Referring  now  to  the  dividing  case  shown  in  figure  lib,  it  is  seen 
that  there  are  several  different  possibilities  for  the  shock  motion,  de- 
pending on  the  initial  condition,  again  indicated  by  circles.  If  the  ini- 
tial condition  lies  outside  the  integral  curves  through  the  saddle  points, 
the  shock  position  merely  oscillates  with  time,  ni'ver  going  through  the 
throat.  If  the  initial  condition  lies  beneath  the  integral  curves  through 
the  saddle  points,  the  shock  moves  upstream,  passes  tlirough  the  throat 
and  disappears;  then  for  the  same  reasons  mentioned  in  the  pri'vious 
case,  it  forms  at  the  throat  at  the  time  corresponding  to  the  first  sad- 
dle point  after  its  disappearance.  It  then  follows  the  integral  curves 
through  the  saddle  points,  so  that  thereafter,  it  just  moves  to  tlu-  throat 
and  never  passes  upstto’am;  (his  motion  is  illustrati'd  in  ligori'  1 5b. 


If  the  initial  condition  should  lie  on  an  integral  curve  through  a saddle 
point,  the  shock  position  is  completely  described  by  integral  curves 
through  the  saddle  points;  the  shock  never  moves  upstream  of  the  throat. 

Finally,  referring  to  figure  11c,  there  are  again  several  possible 
initial  conditions.  If  the  initial  condition  lies  above  the  integral  curve 
entering  the  saddle  point,  the  shock  motion  approaches  a periodic 
form,  never  reaching  the  throat.  If  it  lies  on  an  integral  curve  below 
the  curve  entering  the  saddle  point,  it  moves  upstream  through  the 
throat  and  disappears,  forms  at  the  throat  at  the  time  corresponding 
to  the  first  saddle  point  after  its  disappearance,  and  then  moves  away 
from  the  throat  and  approaches  a periodic  motion,  never  approaching 
the  throat  again.  This  motion  is  shown  in  figure  lie.  Finally,  if  the 
initial  condition  should  lie  on  the  integral  curve  entering  the  saddle, 
the  shock  wave  moves  to  the  throat,  moves  away  immediately  on  the 
integral  curve  leaving  the  saddle  point,  and  approaches  the  same 
periodic  motion  mentioned  above. 

The  numerical  examples  shown  so  far  (e.g.  , figure  12)  have 
been  for  simple  wall  geometries  for  which  it  is  possible  to  derive  an 
approximate  relationship  between  the  parameters  for  the  special 
dividing  case  shown  in  figure  11.  (Eqn.  65).  For  general  geometrii-s, 
it  is  necessary  to  integrate  equation  (4  7)  numerically  along  the  inte- 
gral curves  leaving  the  singularity,  using  equation  (56a)  to  find  starting 

values  near  x =0.  to  find  which  case  holds.  Example's  of  such  cal- 
so 

culations,  for  the  same  wall  shapers  usc'd  in  example  calculations  for 


calculations  made  for  C,  , 0,  that  is.  for  the  case  where  the  steady 

<dd  ’ 

state  solution  is  that  for  which  the  flow  goes  through  sonic  velocity 

at  the  throat  but  is  subsonic  thereafter,  with  no  shock  waves.  Clearly, 

the  unsteady  motion  is  that  illustrated  in  figure  11a.  In  figure  (15), 

two  examples  are  shown  in  which  the  only  parameter  varied  is  the 

steady  state  shock  position,  x^.  Referring  to  the  integral  curves 

through  the  first  saddle  points,  it  is  seen  that  for  x^  1.  5,  the 

situation  is  that  illustrated  in  figure  11c,  while  for  x 0.  75,  it  is 

o 

that  illustrated  in  figure  11a.  Also  shown  in  figure  15  are  the  solution 
curves  from  the  initial  condition  to  the  point  where  the  shock  passes 
through  the  throat.  With  these  two  curves  and  those  leaving  the  first 
saddle  point,  one  can  find  then  the  resulting  shock  wave  motions  cor- 
responding to  figures  13a  (x-  = 0.  75)  and  (lie)  (x  =1.  5). 

U o 

In  summary,  the  above  examples  illustrati'  the  remarkably 
varied  shock  motions  governed  by  the  deceptively  siir>ple  first  order 
nonlinear  equation  (47).  Only  an  outline  of  the  possible  shock  motion 
histories  has  been  presented  here.  In  view  of  possible  applications  to 
inlet  buzz  and  flutter  and  surge  problems  in  turbomachine  ry,  it  appears 
that  more  work  is  called  for  in  this  problem. 

(Ill)  SEPARATION  IN  UNSTEADY  FLOW  FIELDS 

i'he  basic  idea  behind  this  proposed  work  was  to  attempt  to  us»' 
work  done  on  steady  shock  wave  boundary  layer  interactions  at  incipi- 
ent separation  in  ascertaining  under  what  conditions  sucIt  calculations 
could  be  used  in  unsteady  flows.  That  is,  by  ordering  the  partial  tiine 

V) 


m 
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derivalivf  terms  relative  to  the  important  terms  retained  in  the  govern- 
ing equations  in  the  interaction  region,  it  was  felt  that  one  could  derive 
the  conditions  under  which  these  partial  time  derivative  terms  may  be 
neglected.  In  that  event,  the  solutions  found  for  the  steady  flow  could 
be  used  for  the  unsteady  flow  cases  at  each  instant  of  time;  that  is, 
the  boundary  conditions  are  then  time  dependent,  so  the  unsteady  flow- 
solution  is  given  by  a series  of  steady  flow  solutions,  each  with  differ- 
ent external  conditions.  After  finding  these  conditions,  they  were  to 
be  tested  modestly,  by  ascertaining  if  the  parameter  range  in  which 
time  derivatives  are  important  contained  the  case  illustrated  experi- 
mimtally  by  Meier^^\  That  is,  since  Meier's  experiments  with 
shock  wave  induced  separation  in  a nozzle  flow  showed  that  the  separa- 
tion point  motion  lagged  the  shock  wave  motion,  it  was  believed  that 
the  terms  involving  partial  time  derivatives  would  be  impoilant.  It 
has  been  found,  however,  that  except  for  extremely  small  character- 
istics times  associated  with  the  unsteady  flow  external  to  tlie  boundai\ 
layer,  these  time  derivatives  are  not  important.  This  can  be  illus- 
trated by  considering  the  following  two  dimensionless  terms  from  the 
equation  of  motion 

du  , 9u 

a T ^ a X 

where  u — - is  a term  found  everywhere  in  the  interaction  region  except 

a X 

in  the  region  nearest  the  wall.  It  is  of  interest  here  to  comenl  rati'  on 

^^’Weier,  G.  K.  A.  , "Shock  Induced  Flow  Oscillations,  " ACiARn-CF- 
168,  Flow  Separation,  1976,  pp.  .19-1  to  29-9. 
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a partii  ular  rogion,  say  Iho  velocity  defect  layer  in  the  turbulent  boun- 
dary layer.  In  this  r(>gion,  it  can  he  shown*  ^ that  du/dx  = 0{e  ^ ) 

where  f is  the  ordt'r  of  the'  first  time  dependent  term  in  the  velocity 

evaluated  in  front  of  the  shock  wave,  (I.  e.  , u - 1 t c C (t)  + . . . 

u 

Hence,  if  T = t t,  then  for  the  condition  du/  d’l'  « u du/dx. 


T >>  ( 1 69  I 

Clearly,  this  condition  is  met  for  all  but  very  small  values  of  T , 

c h . 

(6 ) 

Moreover,  in  evaluating  the  t value  associated  with  Meier's 

exp»‘ rimeni  s , since  1,  i F>  tnm,  a - 1-10  m/sec,  'I  , ^ 1 msec., 

c h 

-I 

and  r 0.  17,  it  is  si-en  that  t 67  0(k^f  ).  Apparently,  then  the 

I'ange  of  t values  lovered  by  his  experiments  lies  in  the  slowly  varying 
time  regime  and,  indeed,  falls  into  the  regime  referred  to  as  case  |2| 
in  this  report.  Evidently,  then,  (he  time  derivatives  are  not  important 
in  this  cast;;  the  steady  solutions  should  be  valid.  Ihi;  variation  with 
tinu*  of  the  distance  between  the  shock  wave  and  the  separation  point 
evidently  depends  on  the  fact  that  the  shock  motion  is  large  and  hence 
th(?  variation  in  shock  wav'c  strength  is  large  enough  to  cause  this  \ ari- 
ation.  In  any  event,  it  does  not  appear  that  the  cause  is  diu'  to  unsti>ady 
c'ffects  within  the  interaction  ri?gion,  as  originally  suspected.  This  work 
will  be  pursued  to  the  point  of  indicating  whetlu'i-  llie  rather  strong 
timi-  depencience  of  the  distance  tietwei'n  thi'  shock  wave  and  the  si'para- 
tion  point  can  be  predicted  using  known  steady  solutions  for  the 


^Adamson,  Jr,  T.  C.  and  Messiter,  A,  F.  , "Normal  Shock  Wave-Turbulent 
Boundary  Layer  Interactiont  in  Transonic  Flow  Near  Separation", 

Transonic  Flow  Problems  in  Turbomachinery  (Eds.  Adamson,  Jr.,  T.C., 
and  Platzer,  M.  F.  ) Project  SQUID  Report  MiriI-16-PU,  1976,  pp.  392-4H. 


shock  boundary  layer  interaction  and  the  case  (2)  shock  wave  motion 
described  in  section  II  of  this  report. 
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Figure  3.  Lsotac  lis  and  -.■.all  anrl  < enlerliRp  pressure  distributions  for 
for  ronditions  given  in  equations  (9). 
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Figure  3.  (c)  t - 0.78726 
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Fiuurp  4.  r.ent<-rline  pressure  distributions  at  various  times  ffi  r case  (I)  tor  the 
same  conditifjns  used  for  figures  3,  equations  (*^l. 


icure  ( . Isotachs  and  v.all  and  . enterline  pressurr  distributions  for  case  IP.),  for 
- 100  and  all  other  conditions  as  ^iven  in  equations  (9). 
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Fiturp  Comparison  of  case  (1)  solutions  (etjuations  (6)  and  (50b))  and  case  (3),  unified 
■iolutif  ns  (equations  (5Zl)  for  shock  wave  position,  Xg,  for  € 0.05  and  all  other 

■onditions  as  in  equations  (9). 


{•'igiirj'  1 1.  Skf’lrlu's  shovsin^  tiu'  t h m'  possihlo  conti  rat  ion  s lor  intofjjral 
l urvcs  passim;  ttirou^ih  the  saddle  points  (solid  lines  1;  other 
inte^r.il  ^•urves  are  indiealed  h\  dofte<|  lin<-s.  (a)  Inteural  rnr\’es 
leavim:  the  s.iddle  (>oint  rea  i li  tin-  time  axi-^  lietore  ihc'  next 
saddle  point.  (1)1  Inte>>ral  l urves  le.ivinn  the  sadille  point  reai  ti 
the  time  axis  at  tht'  next  saddle  point.  (>  ) Inte^r.il  curves 
leaving  the  sadfile  point  never  return  to  the  titiie  ,ixis. 


Fitiurf'  13.  Sketches  of  shock  wave  motion  v>.hen  the  amfjlilUfle  of  the 

i mp  f(>s  sefl  pressure  o.scillation  is  lar^e  enough  to  drive  th 
shot  k wave  u|>stream  of  the  throat,  for  eat  li  of  tlie-  thret' 
tases  sht)\\  n in  fi^;ure  1 1;  tases  labeled  a,  h,  and  t reler 
to  the  t <)  r re  s pond  in  ^ cases  in  figure  11.  In  each  tase, 
refers  to  the  initial  condition  for  the  shock  position.  (See 
eciuation  15.) 
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Those  curves  which  originate  with  an  x greater  than  any  x on  this 

so  so 

periodic  curve  will  approach  the  periodic  curve  asymptotically  from  above. 
This  periodic  curve  is  nearly  symmetric  about  x^,  the  steady  state  value 
of  X .In  the  dividing  case,  shown  in  figure  11b,  the  curves  entering  and 
leaving  the  saddle  points  are  the  same  curve. 

The  integral  curve  map  obtained  in  any  given  case  depends  upon 
C^d*  1^2’  wall  shape,  f(x).  Although  general  solutions  from 

which  a general  criterion  for  the  dividing  condition  (figure  lib)  could  be 
derived  are  not  available,  an  approximate  result  can  be  found  for  G as 
given  in  equation  (54)  and  f(x)  as  in  equation  (57a).  Then,  equation  (47) 
becomes 


K x 

so  dt 

3 

= -C_-rx  -G  sinbt 

2d  so  n 

(59a) 

^ - T (vTt) 

(59b, c) 

and  the  slopes  of  the  integral  curves  at  the  saddle  points  are  given  by 

equation  (58),  with  t and  G related  as  in  equation  (55).  Now,  if  it  is 

o o 

assumed  that  the  integral  curve  which  passes  through  the  saddle  point  at 

bt  = bt^  and  also  through  the  next  saddle  point  at  bt  = bt^  + 2ir  (e.  g.  , see 

figure  11b)  is  approximately  symmetric  about  bt  = bt  + ir,  then  the  maxi- 

o 

mum  value  of  x^^  is,  from  equation  (59a), 

(X  ) =(-2C  (60) 

so  m 4d 

Next,  if  equation  (59a)  is  integrated  first  over  one  period  (e.g.  , bt  = bt^ 

to  bt  = bt  + 2tr)  and  then  over  a half  period,  then  since  x = 0 at  bt  and 
o so  o 
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^80  ~ ~ °*^®  finds  the  following  relations: 


0 = C,^  + r / X dt 


t = b(t  - t )/Ztt 
o 


(61a. b) 


2 -—  cos  b t 
b o 


(61c) 


where,  in  equation  (6lc),  advantage  has  been  taken  of  the  fact  that  the 
3 

integral  of  over  half  a period  is  half  the  integral  over  a full  period 

because  of  the  symmetry  of  Substituting  for  cosbt  using  equation 

(55),  one  finds  from  equation  (61c)  the  following  relation  for  G . for  the 

o 

special  case  (figure  lib): 

(G^  - G^ = bK(x  )^  /4  (62) 

o 2d  so  m ' ' 

where  K is  given  in  equation  (59b).  Although  this  equation  is  useful  in  setting 
a first  approximation  for  G^»  a more  accurate  result  may  be  found  by  taking 
into  account  the  fact  that  the  integral  curve  in  question  is  not  in  fact  symmet- 
ric, but  is  slightly  asymmetric.  In  this  calculation,  it  is  necessary  to 
employ  an  approximate  form  for  ^ cubic  equation  of  the  following 

form  suffices: 

’'so  ^ ^ 

Now,  at  X = (x  ) , where  dx  / dt  = 0,  t is  defined  as  t , where 

so  so  m so  rn 


t = T + 6 
m 2 


Also,  it  is  assumed  that  6 is  numerically  small  enough  that  terms  involving 
2 

6 may  be  ignored.  Then,  from  equation  (59a)  evaluated  at  x = (x  ) , 

so  so  m 

equation  (61a)  with  equation  (63)  used  in  the  evaluation  of  the  integral. 
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equation  (59a)  integrated  over  one  half  period  (t  = 0 to  t = 1/2)  with 

equation  (63)  used  in  integrating  the  X ^ term,  (x  ) evalviated  usins 

so  so  m * 

equation  (63)  and  equation  (64),  one  can  derive  the  following  relations 

for  (x  ) , 6,  and  finally,  G : 

so  m ' o 


, - bK(x  ) 

(G  - C ) = - 6w  (0.207  r (x  ) - C ) 

o td  4 so  m 2d 


(65c) 


where,  again,  K and  r are  defined  in  equations  (59b,  c). 

Ebcample  calculations  of  the  integral  curves  through  the  saddle 

points,  with  the  sinusoidal  forcing  function  given  in  equation  ^54)  and  with 

parabolic  walls  as  in  equation  (57a),  are  shown  in  figure  12;  the  first 

approximation  to  the  special  value  of  for  case  (b),  calculated  using 

equation  (65c),  must  be  modified  using  trial  and  error.  The  calculations 

were  carried  out  by  numerically  integrating  equation  (59a),  using  equation 

(58)  to  find  an  initial  condition  near  x = 0.  In  the  calculations,  b = 2. 

so  ' 

k = 1,  a = (\  + l)/2  = 1.2,  = 1.5,  = (2  f(x)/(Y  t 1))*^^,  and 

3 

= - 2YC  /3,  where  x is  the  steady  state  value  of  x .In  figure 
2d  uo  o so  ® 

12,  the  letters  a,  b,  and  c refer  to  the  corresponding  cases  shown  in 
figure  11.  In  each  case,  only  the  curves  through  one  saddle  point  are 
shown;  the  repetitive  nature  of  the  curves  at  each  saddle  point  is  not  shown, 
for  clarity.  It  should  be  noted  that  the  value  of  t in  figure  12,  referring 
to  the  location  of  a saddle  point,  is  different  for  each  case.  The  centers. 
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which  also  occur  at  different  values  of  t for  each  case,  are  noted  in  figure 
12.  With  the  parameter  values  given  above,  it  was  found  that  for  the 


special  case  shown  in  figure  11b  the  special  value  for  was,  from  equa- 
tion (62),  (G  ) =4.33  and  from  equation  (65c),  (G  ) = 4.77.  The 

o sp  o sp 

value  which  gives  accurate  results  (figure  12)  is 


(G  ) = 4.968 

o sp 


Thus,  equation  (65c)  is  helpful  in  giving  a relatively  accurate  (4%  error) 

first  guess  for  in  another  case,  with  all  other  parameters  the  same, 

but  with  = 0.75,  it  was  found  that  equation  (65c)  gave  an  estimate  with 

an  error  of  6%.  The  curves  labeled  a and  c in  figure  12  were  calculated 

using  G = 5.5  > (G  ) and  G =4  <(G  ) , repsectively.  In  each  of  these 

o o sp  o o sp 

cases,  curves  entering  and  leaving  the  saddle  point  at  t - t =0  are  shown, 

o 

the  behavior  in  each  case  following  that  sketched  in  the  corresponding  part 
of  figure  11. 

The  solutions  shown  in  figure  12  are  for  very  simple  (parabolic) 
wall  shapes.  There  appears  to  be  no  simple  way  of  predicting 
more  complicated  wall  shapes;  in  general,  it  is  necessary  to  integrate 
numerically  along  an  integral  curve  leaving  a saddle  point  to  see  which 
case  occurs  for  the  given  parameters.  Examples  are  shown  later. 

With  the  mathematical  behavior  of  the  integral  curves  through  saddle 
points  understood,  it  is  possible  to  interpret  the  physical  behavior  of 
the  shock  wave  in  each  case.  Referring  to  figure  11a,  for  any  initial 


